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ABSTRACT 

A design is presented for the experimental investigation of 
the transient growth of an ice layer that forms when superheated 
water flows, under conditions of constant nead, between two parallel 
plates. Both plates are insulated but for a portion of one which is 
convectively cooled to a temperature below the freezing point of 
water. A mathematical formulation is proposed for the case when 
the water flow is vertical. Tne experimental procedure and method 


of testing is outlined, 
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CHAPTER 1 
INTRODUCTION 
1.1, Statement of Objectives 

An important class of heat transfer problems involves solidi- 
fication or melting. Deemed free boundary problems, they are char- 
acterized by a moving liquid-solid interface whose boundary condition 
is nonlinear and whose position is not known 'a vriori' but must be 
determined as part of the solution. For internal flow problems the 
difficulty in obtaining an analytic solution is further compounded, 
Since the liquid phase is directly affected by tie interface motion. 
This results in a coupling of the field equations in both phases unless 
one of the phases can be assumed to be at the fusion temperature, 

The study presented in this thesis is the development of an 
experimental apparatus which will be used to investigate the transient 
solidification of water flowing through a parallel plate test section, 
One of the plates is insulated while the other is cooled to temperatures 
below the freezing point. A system, comprised in »art of this test 
section, has a constant pressure drop applied across it. As freezing 
progresses with time the pressure drop in the test section increases, 
while the remainder of the system experiences a decrease in the pressure 
drop and a corresponding decrease in the flow rate. Conditions under 
which an ice barrier can be grown into the fluid stream need to be 
examined and efforts made to predict the pressure drop required to 
prevent the test section from freezing shut. In this regard an 
experimental apparatus has been designed and constructed so that data 


such as the history of the ice-water interface for various test 
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conditions can be obtained, 


1.2 Review of Revelant Literature 

The first published discussion of a heat conauction problem 
involving freezing or melting was by Stefan [1] . Since then numerous 
mathematical studies nave been undertaken as evidenced by the literature 
review of Muehlbauer and Sunderland 2 | . Due to the difficulties out- 
lined in Section 1.1, few exact solutions for the transient frozen 
layer thickness and temperature distributions have been found. Most of 
these are summarized by Carslaw and Jaeger [3] : 

One area of research which until recent vears has received little 
attention involves solidification or melting on a surface exposed to 
fluid flow. Brush [4] in 1916 discussed ice fornation in water mains, 
He noted that the ice thickness is dependent upon the water temperature, 
flow rate and wall temperature, but presented no taeoretical analysis or 
experimental data. Hirschbere [5] Studied flow through a pipe whose 
wall temperature was maintained below the freezing temperature of the 
liquid, Assuming a constant pressure drop maintained between the inlet 
and outlet of the pipe, a relationship between the perameters for steady-~ 
state conditions and a uniform solid phase thickness was derived, 

Libby and Chen (6 | presented an approximate analysis for "short" 
and "long" times for the crowth of a solid layer depcsited on a cold 
Surface by a flowing gas. The solution was based on the integral 
method as proposed by Goodman 7, 8, 9 | » lLapadula and Mueller [10] 
arrived at a simpler solution to the same problem by employing the 


variational technique developed by Biot by 12 | : 
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Yang (13] Studied the unsteady formation of ice on an instantan- 
eously cooled plate subjected to plane stagnation flow of an isothermal 
liquid, A "short" time solution was derived by asyaptotic series 
expansions, while the "long" time solution based on first-order 
deviations from steady-state conditions was obtained by Laplace trans- 
forms, The mathematical procedure is complicated and application to 
other problems would likely increase the complexity. 

Beaubouef and Chapman fia] used a Similarity analysis employing 
a finite difference technique to give the thickness of a solid deposited 
by a flowing liquid on a cold surface as a function of time and position, 
Zerkle and Sunderland [15] investigated the heat transfer and pressure 
drop effects of steady state liquid solidification at the inner surface 
of a circular pipe, Unlike Hirschberg [5] a uniforn solid-phase shell 
thickness was not assumed, It was argued that the axial velocity would 
retain its parabolic profile throughout the cooling resion,. By assuming 
this form, and by use of suitable nondimensional variables the problem 
was shown to reduce to the classical Graetz problem, Experimental 
observations taken with water flowing through horizontal tubes showed 
that free convection may produce considerable d2viaticn from the 
theoretical results. 

Lee and Zerkle ig| then extended Zerkle and Sunderland's earlier 
work to steady-state liquid solidification in a parallel plate channel. 
By assuming a parabolic velocity distribution in the channel and trans- 
forming the heat conduction equation, a solution was cbtained by separa- 
tion of variables. The resulting dimensionless interface distance was 


obtained by successive substitution into the interface equation. 


ae: 


: 


viii a} Ya nat anereot anets oa eotbate [ef nat 
Teoresisoe? ne ' wal unpade Hols @3 betostdue adelq bafags yf eves 
aster aitenen yas Gi heeled egy notauio: aatt “reoda” & bhuptt 
sebastien ne Seend mattotes enld “pra” act ot hew ,ematonsges 
on saefced gi Fentetao row enotzihno> siate-yhest2 coy) snot isiveb 
of toMsohings gre bataatiaens et avubssevy iestiancemen wT . amet 
<“etesiqnoe a62 vaseani fest! biuay eeisony vaste 
poryolanm sfeylans viral tml: & beew [ay aaeest) wth Tauptuegu 
butlaciah tifoe « to eeandSies siz avin of sephetps? sonweelith stial? a . 
«nots tgeq hae eults To. ng izamt 5.25 Serws hico 4 mp twill, esteat? «yi 
srvneeng, be abana? s6r@ off batagtlcer't feei veel mut tre of eed 
opera sean? edt.te netyeatryitine Htuptl erate ysete Wo cts9¥is gow 
fone srsny-bt ine rvetiw s {2} wvsirgey's afl—ni .eake selurvio o te 
bteaw wWiaoley latue ofS 3008 Deupie ctw 31 .Semeen fon ee copadpte 
GAi dwt es WW AOI eutiasa ae tugteuaing 4ffirre >t steven ed) ntade . 
caldorg ats 2elagiiey lanctermmtbenn efésstya to sau ce ene prec? 2rag * 
(esnwnteogad aldore ctames festeds)> aa! of soubor od mepce a6 | 
bonus extys fasngsbrun Ayes gotwnl? waren tie vader enotievisada 
ait ony nat sates elimeteios aabiw ye mtizevnes sev? tan 
= - «ti fvew, foatiovwers 
TRIAL bem 9698S Lobegree maiz faa] wi drel bos saj 
aiiered a ph ielreotteeiine biuph! s2sze~quasa ot Asau 
ee taslay aticateneg § gators gi . 
— Pe steoagee: mabtouerem seed, mld phere), 
enemy dhitianel sapidahier to nots. 
oriranzque yt bentatde 


Siegel and Savino 17] considered a warm flowing medium solidi- 
fying onto a plate convectively cooled on the underside, Three distinct 
analytic solutions were obtained: (a) an analytical iterative procedure, 
Similar to that developed by Adams [13] » which essentially converged 
to the exact solution on the third iteration, but wiich increased in 
complexity with each new iteration; (b) by employing a second order 
polynomial approximation for the temperature in conjunction with Goodman's 
integral technique 7, a; 2] ; (c) and by applying the integral technique 
and approximate temperature profile to the interface equation, Experi- 
mental values of ice thickness and plate surface temperature throughout 
the transient and steady-state growth period generally agreed favorably 
with the theoretically predicted values, except that under some conditions 
the subcooling of the frozen layer during its formation slowed the growth 
rate substantially, To investigate this effect further Siegel and Savino 

[19] again considered a medium flowing over a plate convectively cooled 
on the underside, but included the heat capacities of both the frozen 
layer and the plate in their analysis, The frozen layer growth was 
obtained by integrating the transient heat conduction equation within 
both the frozen layer and wall. This gave coupled integral equations 
for the temperature distributions in both media and an integral equation 
for the frozen layer growth. The solution to this set of equations was 
then iteratively obtained and agreed well with the experimental obser- 
vations, Savino and Siegel [zo] then applied their analytic iterative 
technique developed earlier to a warm liquid flowing over an isothermal 


plate, It was found possible to carry out a fourth iteration which 
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5 
again essentially converged to the exact solution and which showed good 
agreement with Libby and Chen [6 | » Lapadula and !fueller io | » and 
Beaubouef and Chapman ia] ‘ 

Gort 1] studied analytically and experimentally the freezing 
problem inside a vertical cylinder both with and without forced flow. 
An approximate solution was derived for a fluid whose temperature was 
at its freezing point and which flowed through a cylinder whose wall 
temperature was maintained below freezing, 

DesRuisseaux and Zerkle [22| developed a technique which predicted 
conditions under which a hydraulic system freezes shut. The system 
pressure drop was analytically determined as a function of Reynold's 
number at steady-state conditions for various wall temperatures. The 
minimum pressure drop to stop the system from fr2ezing shut was specified, 

Ozisik and Mulligan fe3| investigated analytically the transient 
freezing of a liquid flowing inside a circular tube, The restricting 
assumptions were: 

a) the tube wall temperature was lower tian freezing 

b) slug flow velocity profile 

c) the variation of the thickness of the solid layer with both 

time and position was small so that a quasisteady state heat 
conduction in the solid phase could be assumed, 
The method of solution involved determination of the dimensionless 
liquid and solid phase temneratures in terms of the unknown dimension- 
less solid-liquid interface radius, The equations were coupled to the 
interface equation to obtain a single ordinary differential equation 


for the interface radius which was then solved to give the interface 
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radius as a function of time and axial position, 

Stephan [24] considered flow along a convectively cooled plate 
and through a convectively cooled pipe. A parabolic temperature 
distribution was assumed for the solid phase matcning all the wall 
boundary conditions but satisfying tne energy equation at the solid- 
liquid interface only. An interesting outcome of the analysis was 
that in pipe flow, growth of tie solid layer could stop at two critical 
points for certain values of the ambient heat transfer and the convective 
heat flux at the solid-liquid interface. One of tiese points was a 
Stable condition, the other unstable; both however, were dependent on 
the preceeding history of the solid phase as to their SLavilacy. 

Freeborn [25, 26 presented a theoretical and experimental study 
of the transient growth of an ice layer formed as Superheated water 
flowed through a convectively cooled vertical cylinder, The analysis 
considered two regions: an ice free zone where the fluid's sSuperheat 
prevents ice formation, and an ice forming zone, Following the work 
of Sellars et al [27] the temperature profile at the entrance to the 
ice forming zone was obtained. By assuming the parabolic velocity 
profile used by Zerkle and Sunderland hs | which adjusts for local 
variations in the ice thickness, Freeborn derived an asymptotic 
Solution for the temperature distribution in the ice and was able to 
predict the rate of solidification and its effect on the pressure drop. 

Nyren [28, 29] extended tne work of Freeborn [26] to include the 
effect produced by circumferential variations in the external heat 


flux aS well as the sensible heat in the ice layer, A perturbation 
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7 
ez~pansion about the quasi-steady-state solution was found. The influence 
of the non-uniform convection on the interface profile was shown graphi- 
cally for two cases, Experimental results for flow through a unifornly 
cooled and inclined pipe were given, 

“ore recently, due to the complexity of the problems being investi- 
gated, solutions have tended to rely on numerical analysis. Bilenas and 
Jiji [30 . 31 | presented both a variational and a numerical solution to 
the problem of axisymmetric fluid flow in tubes with surface solidi- 
fication, The numerical solution was based on finite difference approxi- 
mations of the continuity, momentum and energy equations with Lagrangian 
interpolation equations employed at the moving interface. The paper was 
unique since prior to its publication no numerical scheme was available 
which would solve two-dimensional fluid dynamics problems involving 
solidification, Effects such as viscous dissipation, axial viscous 
Shear, axial conduction and discontinuity in the mass density due to 
phase change were included in the analysis. 

Although not concerned with the field equations governing the 
liquid phase, the publications of Lazaridis [32] and Bruno [33] should 
be mentioned since both were concerned with transient multi-dimensional 
change of phase. Lazaridis devised one-, two- and three-dimensional 
finite difference schemes to solve the heat conduction equation for 
typical boundary conditions. Bruno treated the one- and two-dimen- 
Sional heat conduction equations describing a finite body bounded by 


two insulated surfaces subjected to a known heat flux. 
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1.3 Applications 
With the formation of ice and snow and the annual coming of spring, 
solidification and melting were known to earliest man. Today free 
boundary problems are encountered in such diverse areas as the 
freezing of foods, casting of metals, spacecraft propulsion systems, 
and the control of frost formation. Of particular concern to Canadians 
is the industrialization of the country's northern and arctic regions, 
Problems such as building pipelines through areas of permafrost, and 
failsafe heat and water supply systems for arctic cities require a 
thorough understanding of the freezing and melting process, This is of 
special importance if the environment of the region is to be preserved, 
The application that motivated the present study is the possibility 
of controlling and stopping the flow of water by an essentially static 
means. Such a device which could function as a valve in part of a water 
piping system would have the advantage of involving a minimum of moving 
parts, This also reduces the sealing problems associated with the 


penetration of fluid boundaries by shafts or other devices, 
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CHAPTER II 
MATHEMATICAL FORMULATION 

2.1 Governing Equations 

Consider a system where initially superheated water flows vertically 
upwards between two insulated parallel plates with a steady laiinar 
motion, At some time, say t = U, a portion of one of the plates has 
its temperature lowered to a subfusion value. Almost instantly a thin 
layer of ice forms on the surface of the cooled wall producing regions 
referred to as the adiabatic and freezing zones in figure 2.1. As time 
progresses the ice grows towards the opposing side until either bDlock- 
age occurs or a steady-state profile is reached, Uuring the freezing 
process the pressure drop across the system is maintained constant and 
the temperature distribution at the inlet to the freezing zone invariant, 

If the flow through the freezing zone is asSsuined laminar, and the 
liquid Newtonian with constant properties (p, Kis Cis u); the differen- 
tial equations governing the liquid phase in two diilensions are: 
CONTINUITY EQUATION 
“ + = = 0 2eig = 1 
EQUATIONS OF MOTION 
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SOLID PHASE 


FULLY DEVELOPED 
HYDRODYNAMIC 
REGION 
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LEADING EDGE OF 
FREEZING ZONE 


ADIABATIC ISOTHERMAL ENTRY 


DIRECTION OFSELOW 


FIGURE 2.1 PARALLEL PLATE SYSTEM WITH 
SURFACE SOLIDIFICATION 
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If the material properties (p., k,, C, ) for the solid are also 
constant but not necessarily the same as tne liquid phase the energy 


equation for the solid becomes: 


Oise. 9°Ts + 07Ts 


2,c= Order oft Magnitude Analysis 

The conditions under which the foregoing equations may be simpli- 
fied can be obtained by non-dimensionalizing the equations in terms of 
dimensionless variables of unit order. This normalization procedure 
in addition produces dimensionless parameters essential for experimental 
Simulation of the phenomena to be studied, for determination of the range 
of experimental data [34] : 

The first step in the procedure is to normalize the governing 
equations by determining suitable scales or characteristic quantities 
whicn reduce the variables and their derivatives to non-dimensional 
form of unit order of magnitude, These characteristic quantities, chosen 
in accordance with the physical system the governing equations describe, 
form upon nomalization dimensionless coefficients of the terms having 
unit order, Tnese coefficients, or dimensionless parameters then, 
give an indication as to the importance of each of the terms in the 


governing equations, 
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The variables in normalized form become: 12 


t° = t/te - time coordinate 
X7 = xX/X_ - axial space coordinate 
y~ = Y/Yc - transverse (normal) space coordinate 
Owen tL == Pe - pressure 
Pe 
u~ = u/Uc - axial velocity 
Vitec " - transverse velocity 
6, “= us = : - liquid temperature 
6s “= io = : - solid temperature 


Where: p, =the fluid pressure at the freezing section outlet 


{ 


To = inlet fluid temperature 
Tt = fusion temperature of the liquid 
Tw = the wall temperature of the cooled plate 


Tne superscript "“", designating a dimensioniess variable of unit order 
Of magnitude, while the subscript "c", indicates a characteristic 
quantity which is to be determined, 

Using these quantities the continuity equation 2.1 - 1 becomes 


oUG BEM EXCVG ER OVe = 0 


Since the axial and transverse velocity gradients are of the same 


magnitude then 
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Substitution of the variables into the momentum and energy equations 
(2.1 = 2 through 2.1 - 5) leads to the following equations. 


Momentum Equation = x component 
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Momentum Equation - y component 
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Energy Equation - liquid region 
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The governing equations are written in the above form as it is 
expected that in the liquid region the convection and advection terms 
Will be important, and hence of order unity. For the energy equation 
in the solid the normal temperature gradient is a dominant term if 
growth of a solid phase is to occur. 

A suitable characteristic quantity for ye in the solid region 
would be the steady state (fully developed) solid thickness. If this 
value is teemed s,, the time scale t. may be determined from the 
energy equation of the solid by noting thet the transient is of the 


Same order as the transverse conduction only if: 


osCs | se* | x] 
Ks tc 
2 2 
Hence t. = sls mary 7k 22. ep 
2 Os 


The time scale could have been determined from tne governing equations 
for the liquid region. It, however, would have involved determination 
of additional characteristic quantities. It should be noted that the 
above definition gives the rate of solidification of the solid phase 
directly in terms of the properties of the solid, 

Examination of the foregoing equations reveals that after some 


manipulation the following dimensionless groupings may be made to 


appear: 
Cee ea Ole » the thermal diffusivity ratio 
As 
Reg = UcYc » the characteristic Reynold's number 
U 
or ee » the Prandtl number 


aottiupe ene at 
* wrod nae va. p 2 dae 


ae tget pifog Sit nt ov 1c’ yt rea ‘srjetpoatens 
x brfoe (he jotavah aftuh) sfate qpeste eid st plove : 
mit? St, ¢ Someey et ater 


om w 


r ear? pantaneds od, yen —3 
+ gain pifaon Qo bhfae s aly Yo fnoltapps arene 


ey) a 25 ~sbt0 ome 


si Yo et saat Sit? 


enor ait arenes ont iKAT 1 
“notgunt maga havfovet evan biiow ~Wevenon . 


od sane bosin so \iwodd 23 . Sati Soap She 


sent rte Ne notanotint ve 2 ta -gSar 409 covhy nor tattafh ave 


7 
paoe 49 Sine z 


aint 2 zntauore resteohensa 


_ 


oe 


as 


bean Note 


a 
‘J 
7 


st? Vint uotsaubnoo 


ing ‘Aotsaavoos wt? notged Ptuptl al at tats 


_ » Ds 
3) B26 


5: ovav 


1 
e *, pts = 


a 


7 


. 


S407 ,ettne Wobt0 So Bane bees _dasrroml sd 
Fhe aed evgories lesen wit btfos att at 
asd as zt seneig pitoe i to —e 


stiettur A 


i 
a hy, i= ex 
| ; — 
L J 
2 oe = ms 
‘= 
e 


Pee bo 


tniostEMs 


bese re ng mine yusoca et aig 1 psn ema 
i gateaitat ais mobbafuc 


nga 


- _ 


= wens -) 4° 


~*~ 


basooqaa. 


itis 
. 


poet 
gariesiah dase avs pives stase acsks ont 


_ 
if ah) WR 
» Fong hatha w 


aston aie 10 esis canoes oud Oo eave? at qiszent 


7 


a 


‘ 


Peas eRe ee Py » the Péclét number 
Uc 2 

Br? 2s Veg » the Froude number 

Fu = ome » the Euler number 
0, Uc 

oa: Ws » scale factor liquid region 
Cc 

Gg = es » scale factor solid regions 


The governing equations become: 


Momentum equation = x component 
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Momentum equation - y component 
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Energy equation - liquid phase 
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is the dissipation function, 


Energy Equation - solid phase 
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Axial conduction in the solid phase in comparison with conduction 
in the transverse direction is negligible provided 65 <<1. Similarly 
in the liquid phase axial conduction can be ignored in relation to the 
transverse conduction term if 6 << 1. Comparison between the advection 


term and the dissipation term shows that for_E_ << 6, viscous dissipa- 
Re 


tion can be neglected, The value of 6 used, moreover, depends on the: 
choice of the characteristic quantities yand x,. Like the classical 
boundary layer problems this choice depends upon tae region of interest. 
ear the entrance of the freezing zone x —»O so that a suitable 
value for y, mignt be the boundary layer thickness, while for regions 
farther downstream tie choice could be the fully developed fluid 
thickness, 

In a fashion similar to the energy equations, examinations of tne 
zu" , 


xX = component of momentum reveals that the viscous force term 9 
oy2 


is small when compared with the viscous force term represented by  92u’ 
dy2 
provided 6 << 1. Moreover, if Re, >> 1 >> 6, the magnitude of both 


viscous forces would be negligible, though the term 32u’ would still 


oy? 
have to be retained in order to satisy the transverse boundary 


conditions, As the flow is assumed vertical the only body force which 
occurs is in the axial direction, It may be disregarded provided 


_l<< 6. Experimentally, with the low velocities which would be 
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incurred, this condition did not prevail. The driving force, 0, 
OX” 


was naturally also retained, implying Eu = 1 from which the scale factor 
for the pressure is found to be : Po PUc Tne y = component of 


the momentum equation reduces to dp~= 0 provided << 1. 


a 
oy” Ree 

2.3 Mathematical iodel 
Based on the results of the previous section it is advantageous to 


introduce the following set of dimmensionless variables. 


= a t ° fe ° . 
“§ “Se dimensionless time coordinate 
ao & al dimensionless thickness of the liquid 
X = x dimensionless axial coordinate 
u 
Yas y dimensionless transverse coordinate 
P = Deas dimensionless pressure 
olm2 
U = 7 dimensionless axial velocity 
m 
Vo aan dimensionless transverse velocity 
6 = USS: dimensionless liquid phase temperature 
S Rani 
6. = eee wih dimensionless solid phase temperature 
S Vast 
Where d = the distance between the parallel plates 


Um = the average velocity at time t = 0, 
The mathematical formulation defined by the simplified governing 


equations, initial and boundary conditions becomes: 
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Momentum Equation 


X = component 


OOP, OT d oY On 6Re OY? OFy 
Y = component 
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Energy Equation 
Liquid Phase 
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With the above initial conditions and the following conditions, 
U(0, 0.7) = 0 


U(O, 1,7)e= 0 
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Wall Boundary Conditions 
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Interface Boundary Conditions 
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Examination of the above reveals the following dimensionless 


parameters, 

6 = a scale factor 

a = OL thermal diffusivity ratio 
As 

Re = Und Reynolds number 
Wy) 

Pe = Umd Péclét number 
o7y 

= LS ratio of thermal conductivities 
Ks 

o = Te - Tw Superneat ratio 
To = Tf 


Ste= Calle salt) Stefan number; ratio of liquid sensible heat 


to the latent heat 


Fr = Um? Froude number 
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In order to ensure that the velocity profiles obtained from the 
momentum equation 2,3 - 2 are valid, the pressure gradient term must 


Satisfy the integral continuity equation. 


1 n X 
UdY = { var = l-o | an dx 2.3 <9 
d5aP oT 
0 0 9) 


The equation is derived in Appendex A by considering a control volume 
Surrounding part of the freezing section. The first term representing 
the axial mass influx at X = 0, while the second term is the axial inass 
efflux at some distance X. The difference in these terms is the mass 


flow at the liquid - solid interface, 


2.4 Assumptions 
The assumptions used in the derivation of the mathematical model 
represented in section 2.3 may be summarized as: 
1, The liquid flow is assumed to remain laminar and two 
dimensional throughout the freezing section. 
2. The liquid is assumed to be Newtonian possessing constant 
viscosity. 
3. Both the solid and liquid phases have constant density, 
constant thermal conductivity, and constant diffusivity. 
4, The temperature of the cooled wall is assumed constant and 
below the freezing point of the liquid, whereas the Opposing 
wall is assumed to be adiabatic. 


9. Axial conduction in both phases is negligible, 
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In the liquid phase the viscous force due to the axial 

Shear, dissipation and transverse pressure gradient are neglected, 
The velocity profile of the liquid at the entrance to the 
freezing section is fully developed but not invariant. 

The liquid temperature at the inlet of the freezing section 

is constant, 


The pressure drop across the system is constant, 
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CHAPTER III 
EXPERIMENTAL INVESTIGATION 
3.1 Requirements 
It was intended that the experimental data obtained be representative 
enough so that a comparison with any analytical study of the model pro- 
posed in Chapter II could be made, The apparatus, moreover, was to be 
sufficiently flexible that other conditions could be investigated, in 
particular the effect of turbulence and natural convection. The chief 
concern of any study involving a moving liquid undergoing a solidification 
Or melting process is the amount and rate at which the phase change occurs 
and the effect on the system pressure drop. This imposes the obvious 
difficulty of determining the location of the liquid-solid interface 
without disturbing the flow. 
The experimental design as a physical realization of a mathematical 
model can only hope to approach the conditions the model prescribes. 
It must, nevertheless, accurately produce the assumptions imposed by 
the theoretical analysis if a comparison is to have validity. This 
implies that several constraints be applied to the experimental design 
and testing procedure: 
1, Since the test section must have a finite width, any channel 
chosen would have an aspect ratio, Y, which approached zero. 
2. The hydrodynamic entrance length was to be sufficient to 
produce a fully developed velocity profile at the inlet to 
the freezing zone for all laminar flow conditions, 
3. Boundary layer formation on the side walls was to be negligible, 
4. The surface roughness, waviness, discontinuities, and section 


variances were to be minimized as best as possible, 
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5. Axial conduction resulting from the step temperature change 
was to be minimized, 
6. The growth of the solidified layer was to be of sufficient 
thickness to obtain accurate results, 
7. The material for the insulated walls was to be a poor enough 
heat conductor to eliminate ambient heat conduction effects, 
8, None of the dimensionless parameters should violate any of the 
assumptions posed by the order of magnitude analysis. 
9, The total system pressure drop should not vary appreciably 
during a test, 
In addition, it was felt desirable to observe the solid layer during 
the transient conditions. This necessitated the side walls, at least, 
to be constructed of an optically transparent material. Although some 
of the above requirements apnear to conflict, an adequate design can 


still be properly constructed, 


3.2 Description of Exnerimental Apparatus 

The apparatus was considered to consist of two separate systems; 
the "coolant" system which determined the boundary condition on the cooled 
wall of the freezing zone and the "warm water" system wnich provided the 
internal flow between plates. A schematic of the apparatus is snown in 
figure 3.1, photographs of the system in figure 3.2 and 3.3, while con- 
struction drawings are given in Appendix B, A list of manufacturers of 
the various equipment used is given in Appendix C, 

The coolant system was a closed aa consisting primarily of a 
reservoir, pump, control valve, flexible piping and test section 


cooling jacket. The coolant fluid, a water-etiylene glycol mixture, 
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27: 
was pumped through the cooling jacket where it gained heat transferred 
from the “warm water" side of the freezing section, The coolant's 
temperature was maintained by means of a temperature controlled three 
way valve which mixed recirculated coolant with refrigerated coolant 
from the reservoir, 

The "warm water" system was an open loop where water at a constant 
temperature was pumped from a temperature controlled storage tank tc an 
overflow reservoir, Lxcess water from the reservoir spilled back to the 
Storage tank while the remaining water flowed to tie test section. This 
water was then forced upwards between two insulating parallel plates 
allowing the velocity profile to become fully developed before entering 
the freezing zone, As it passed through the freezing section the water 
cooled and some of it solidified forming an ice layer on the cooled wall. 
The cooled water continued to the upper reservoir where it spilled over 
a weir and flowed down into a turbine flowmeter from which it finally was 
discharged into a drain. Since gravity acted as the driving force, the 
total head acting on the system was denendent on the height differential 
between the top of the weirs in the overflow and upper reservoirs. Thus 
the initial flow rate was varied by raising or lowering the overflow 


reservoir, 


Oren lve est oeECt1on 

The test section or freezing zone consisted primarily of two 
rectangular channels separated by a heat conducting wall. A simplified 
cross sectional view is shown in figure 3.4, A photograph of the "warm 
water" side may be seen in figure 3.5 while the "coolant" side is shown 


in figure 3.6, For a more detailed description tne reader should refer 
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to Appendix B, DWG Bl to DWG B4, 

The dividing wall was % inch thick copper plate which provided a 
freezing length, 1, of 24 inches and a freezing wide, b, of 10 inches, 
Copper was chosen because of its high thermal conductivity, low heat 
Capacity. The plate was clamped between channels constructed from 1 
inch and 1; inch thick plexiglass, Plexiglass was chosen as it provided 
good visibility and insulation properties, Leakage between the "warm 
water" and "coolant" sides was prevented by use of 1/8 inch O-rings, 
and by chloroform welding of the plexiglass channel joints. 

The distance, d, between the plexiglass and the copper plate on 
the "warm water" side was set at % inch, which gave an aspect ratio, 

Y, of 0.05, Rudimentary calculations indicated that with this aspect 
ratio, the effect of boundary layer crowth on the side walls could be 
ignored for all but the lowest Reynold's numbers, 

To ensure efficient heat transfer from the coolant, flow in the 
coolant channel should be turbulent, Due, however, to the high viscosity 
of the coolant at temperatures in the operating range - 30°F to 10°F, 
and also due to a limited pump capacity; the coolant flow could not be 
made turbulent without resorting to very small gap sizes. Instead the 
depth of the coolant channel was made 1 inch to suit space limitations 
required by the thermocouple probes. Provision was made for vortex 
generators in the form of triangular plexiglass blocks to be cemented in 
a staggered pattern to the plexiglass wall. Since the inlet flow to the 
coolant channel was turbulent (because of the small diameter of the inlet 


piping), the generators acted merely to sustain the turbulent flow, 
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3.20ee eadiabatic Section 


Tne section was in essence a continuation of the "warm water" side 
of the test section except that it was completely constructed of plexi- 
glass, The upper 13 inch thick plexiglass plate was continuous the 
length of the apparatus, while the lower plate was shorter butting 
against the copper plate of the freezing onerienl A smootn transition 
was made from the plexiglass to the copper plate. 

The adiabatic section provided the entry length required to ensure 
the velocity profile entering the freezing section was fully developed. 
From the study done by Han [35] the entry length required was calculated 
to be 28,6 inches, The section, however, was constructed with a 32 
inch length. Dye injected into the fluid stream as a check showed that 
the velocity profiles were in fact fully developed at the entrance to 
the freezing zone. 

3.2.3. Lower and Upper Reservoirs 

The purpose of these reservoirs was to mainly provide means of 
supplying and discharging fluid from the test section. The lower 
reservoir acted as a diffusing section for the fluid entering the 
adiabatic section, This was useful as it reduced the head loss that 
would have been caused by a sudden area cross section enlargement. The 
upper reservoir served as the discharge point for both the "warm water" 
and the "coolant" system. A weir 8 inches in heignt was used as the 
discharge point for the "warm water" system, 

The reservoirs were entirely constructed of plexiglass except for 
the piping connections, Joints were welded by chloroform as in other 
portions of the apparatus; 1/8 inch O-rings provided the necessary seals. 


The construction drawing of the upper reservoir is given in DWG = B5 
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while the lower reservoir in DWG = B6, both in Appendix B, 


3.2.4, Overflow Reservoir 

The overflow reservoir consisted of a rectangular plexiglass tank 
mounted on a carriage constructed from 1 inch angle iron, The carriage 
utilized bal] bearing wheels which were constrained to move vertically 
by two 7 foot 6 inch channels bent from 16 gauge galvenized sheet metal. 
A pulley system consisting of a pair of pulleys, windlass and steel wire 
Strand was used to raise or lower the reservoir and carriage. 

The reservoir itself was composed of two compartnents separated by 
a weir 14 inches in height. "Warm water" was pumped rapidly enough from 
the storage tank to fill the first compartnent, Some of this incoming 
water was discharged from a pipe connection at the base of the compart- 
ment and flowed towards the test Section, The rest of the water Spilled 
over the weir into the second compartment where upon it was returned to 


the storage tank, 


3.2.5, Apparatus Support Stand 

The adiabatic and freezing sections were Supported by a stand 
constructed from 3 inch by 2 inch rectanaular tubing, The sections were 
bolted to a mounting plate which in turn was bolted to an axle. The 
axle was Supported by two roller bearings mounted onto two arms of the 
Stand, A gearbox keyed onto the axle shaft allowed the test section to 
be rotated to any angle 0°to 90° off the ver ticalos aiWils: was accomplished 
by turning a hand crank on the gearbox, A counterweight system was used 
to reduce the loading on the worm-gear drive, Although the gear ratio 
was high enough so that no additional braking was necessary, a locking 


mecnanism was added to eliminate chance movement, A simple needle 
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pointer arm, Screwed to the centreline of the axle and extending over 
the protracted range, indicated the amount of inclination of the 
apparatus, 

Prior to testing the stand was levelled and bolted to the 
laboratory floor. The gearbox, pointer and locking mechanism can be 
seen in figure 3.7, A photograph of the apparatus shown inclined is 


given in figure 3,8, 


3.2.6, Instrumentation 
The primary measurements recorded during the course of an experi- 
mental test would be: 
1) the ice layer thickness throughout its transient growth and at 
its steady-state position 
2) the surface temperature along the copper plate as a function 
of time 
3) the bulk temperatures of both the coolant and water flowing 
the test saclay 
4) the mass flow of the water as a function of time 
5) the pressure profile tiroughout the test section as a function 
of time 
While most of the measurement techniques used were standard and call for 
little explanation, others may require a more detailed description, 
Perhaps the most difficult instrumentation problem encountered with 
solidification studies involves the measurement of the solid phase thick- 
ness, Siegel and Savino [17] employed three techniques referred to as 


photogranhic, visual (utilizing a cathetometer and grid Sheets), and 


~tiaqas gm Wee wi crib Sabraney etmaamuesonn etambas aT 
a sy tat lawom, 2989 fadeae 
£ bas adr teatames: crt benindereneNet Nk aciad e 

- wnfatang axetevbest a Rat | 

ncttaqm® § @e agete woses <n asnenncilanen tt - 
Z ante Ye -% 

onchiaina eas instoos std ot Ye codecs dda tb 

seats 


rn 


votsoeut & ap ABERDbs teed ons saonnuaya al the azer ee 


_ - a ow ona 


ae 
a ; ‘ i 7 
2 ow se oui Rapti 12 to z 
7 ve : : 
7 7 - yy, 
mh i 4 7 , > 


34 


NOIL99S 1841 G4ANITDNI 8°€ FIND I4 


YILNIOd 9 XOGYVID Z 


€ AWdNDI4 


ye 


i 
i 
: 


. 12 


. j 
11 ; 7 
€ 
ad 


35 
temperature probe methods, Gort [21] developed a probe which deployed 


mechanical arms to detect the solid phase Surface, Nyren [29] used a 
finger like sensing arm and voltage transducer to produce a continuous 
interface profile reading. Bailey and Dula [36 ] and Freeborn [26] 
derionstrated that ultrasonic signals can be made to accurately determine 
the water-ice interface without the inherent disadvantage, associated 
with mechanical probes, of disturbing the liquid phase, 

Ultrasonics were adopted as the method seemed to provide the most 
feasible solution, Due to the narrow gap between the test section wall, 
any mechanical probe devised would have required a high degree of 
miniaturization. The problem, moreover, seemed ideally suited to the 
acoustical technique as copper has an exceptionally hich reflectivity to 
sound waves, This enabled a clear distinction to de made between sound 
Signals returning from the various interfaces, 

The ultrasonic unit, a Branson "Sonoray 600", was operated in the 
pulse-echo mode using a single transmission-recciver transducer, The 
transducer probe would generate an initial pulse (usually referred to as 
the big bang). Reflections, or echos, would be received from each sur- 
face as shown by the typical pulse-echo Signal depicted in figure 3,9, 
The time lag between the echos returned from tne copper plate and the 
ice-water interface determined the thickness of the ice, Through the 
use of an electric timing gate, the Branson Tine Analog Gate, an analog 
Signal of the ice thickness was received, This Signal corresponded to 
a precalibrated thickness, 

The probe was held flat against the outer plexiglass surface of the 
test section by a carriage, as illustrated in figure 3,10, A lead screw 


connected to a reversing D, C, motor allowed the carriage to traverse 
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Vee PULSE (BANG) FROM TRANSDUCER 


ECHO FROM PLEXIGLASS WALL 
ECHO FROM ICE-WATER INTERFACE 
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FIGURE 3.9 TYPICAL ULTRASONIC PULSE-ECHO 
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back and forth along the centreline of the test section. A thin film 


of oi] was needed to act as a couplant between tie transducer and 
plexiglass surfaces, A ten turn potentiometer monitored the position 
of the probe by varying the voltage from a fixed datum, Signals from 
the timing analog and the potentiometer were then used to drive a X-Y 
recorder, Thus a profile of the ice-water interface throughout any 
portion along the centreline of the test section could be generated, 
In addition, the time required to complete a traverse was variable, 
and depended either on the power given to tne D. C. motor or the 
positions chosen for tne limit switches, This proved to be a useful 
feature as higher solidification rates necessitated faster profile 
generation in order to record the transient growth, A photograph of the 
ultrasonic unit used is given in figure 3.11 .— while the traversing 
unit is photographed in figure 3,12, 

The surface temperatures of the copper plate were measured by a 
series of specially constructed thermocouple probes snaced throughout 
the length of the plate, The arrangement of tie probes concentrated 
most of them along the centreline of the plate and near the leading 
edge where most of the transient effects occurred: see DWG - Bly 
Appendix B, A typical probe (DWG - B2) consisted of a pair of 36 gage 
copper-constantan wires soldered to a copper insert disk. A surface 
groove 0,010 inches deep, 0.010 inches wide and 0.125 inches in length 
was tmachined into the insert and the thermocousle wire passed through a 
1/16 inches diameter stainless steel tubing and soldered to 28 gage 
thermocouple extension wire. 

The purpose of the thermocouple probe design was to provide accurate 


Surface temperature measurements, particularly during the transient 
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conditions, The small diameter Provided very rapid response times with 
a minimum of surface distortion, Initially attempts were made to fix 
the thermocouple wire directly to the copper plate, Due to the high 
heat conductivity of copper and the large surface area involved, 
Soldering was not attempted, It was feared the amount of heat required 
would have induced excessive thermal distortion of the plate, Instead, 
Surface grooves were machined into the plate and the thermocouple wires 
with their tips bared of insulation laid in place. The Surrounding 
Copper was then peened over the tips providing good thermocouple 
junctions, The wires were passed through holes drilled in the plate . 
and the holes and remaining portion of the surface grooves filled with 
epoxy. Leakage, however, from the "coolant" to “warm water" side 
developed after a few preliminary runs. Due to the thermal expansion 
and the high fluid pressure on the "coolant" Side, the epoxy was unable 
to adhere to the copper surface, 

A particular advantage gained by the use of tie redesigned 
thermocouples is that each thermocouple was a removable probe 
individually calibrated prior to installation (Chaster Iv), 

To determine the bulk temperatures of the fluids a standard 
sesitens toh type probe was constructed, The probes were located at the 
inlets and outlets of the test Section, and as well in the adiabatic 
section, Locations of the probes are given in Appendix B, DWG - B61, 

The “warm water" system temperature was set and maintained by two 
thermostatically controlled 1500 watt immersion heaters installed in the 
constant temperature storage tank, or by a & ton refrigeration unit 
(figure 3.1), The temperature of the "coolant" System, however, was 


controlled by a resistance bulb thermometer located in the inlet line 
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4] 
to the test section, The thermometer produced a signal proportional to 


the desired temperature which activated a motor operated three-way control 
valve, This allowed previously cycled coolant to mix with the colder 
coolant from the reservoir, Specifications for both systems has been 
given in Appendix C, 

The mass flow of the water flowing through the test section was 
obtained by use of a turbine flowmeter, The meter was located after 
the test section as there the pressure drop would not have any effect 
on the upstream conditions, A pickup generated an output frequency 
Proportional to the flow rate. This sianal was converted to a dince = 
milliampere output used to drive a continuous chart recorder. As the 
meter was accurate only between the range 0.8 = 8 g.p.m. an alternate 
collection system was used (by by-passing the turbine meter) if 
necessitated by the flow rate, 

The pressure profile throuchout the adiabatic and freezing zones 
Was obtained by a series of static holes drilled along the outer 
plexiglass wall, The leads from these holes were connected to a multi- 
tube inclined mercury manometer. Ouring the course of an experimental 
run a 35 - millimeter single reflex camera coupled with an automatic 
timer could photograph the manometer at predetermined time intervals to 
record the transient pressure variations, In addition to obtaining the 
pressure profile, two static holes were placed in line with the tops of 
the weirs in the upper and overflow reservoirs, Leads from these »ressure 
taps were connected to a nigh sensitivity capacitance differential 


pressure transmitter, 
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CHAPTER IV 
EXPERIMENTAL PROCEDURE 
4.1 Calibration 

Preliminary to actual testing the various instrumentation used by 
the apparatus required calibration, The Honeywell Electronik 15, low 
sensitivity twenty-four channel recorder was calibrated in conjuntion 
with the thermocouple probes prior to installation in the test section, 
The probes were immersed in a Rosemount Engineering Model 910AC cali- 
bration bath whose temperature was monitored by a Hewlett Packard quartz 
thermometer model 2801A, which had been previously calibrated against 
a platinum resistance thermometer, A calibration curve for each of the 
probes was determined for temperatures within the rance of - 30°F to 
100° FS) the factory calibration of the Brooks model 4901 - 10 turbine 
flowmeter was checked by collection of the water discharge, The 
Rosemount Engineering iiodel 1151 DP differential pressure transmitter 
was calibration against an inclined water ianometer. 

In order to obtain meaningful results from the ultrasonic thickness 
tester (Branson Instruments Sonoray 600) a means of calibration was 
developed. A small model duplicating the test section configuration 
Was constructed and ice of a uniform thickness grown, The thickness of 
this "calibration" ice was measured by means of a micrometer probe 
similar to that developed by Sieqel and Savino [27 ] . Essentially it 
consisted of a depth micrometer from which a taermocouple wire probe 
extended, The junction of the thermocouple was carefully lowered to 
the vicinity of the ice-water interface while the temperature response 
Was noted, The position of the interface was considered determined 
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when the temperature reached the value of the freezing point. Using 

the ultrasonic tester the thickness was measured and adjusted to agree 
with the value obtained from the micrometer probe. The ultrasonic unit 
Was now considered to be calibrated; various steel gauge blocks were 
measured and the results recorded, Thus as part of the experimental 
test procedure, the ultrasonic unit's calibration would be readily 
checked by again measuring the gauge blocks, and if necessary adjusting 


to reproduce tne previously recorded results. 


4.2 Test Procedure 

The tests wnich can be conducted can be thought to be of two types, 
no-flow tests and flow tests. \No-flow tests are performed to indicate 
the reproducibility of the results and to provide a cohiparison witn 
Studies available in literature. The procedure to be followed in con- 
ducting a flow test (a no-flow test having obvious deletions) would be 
as follows: 

1, Water would be added to the constant temperature storage tank, 
and the temperature control set as desired. 

2. The refrigeration unit would be turned on to cool down the 
coolant reservoir. 

3. After a prolonged period, sufficient to allow for the escape of 
dissolved gases, and long enough for the water temperature to 
have stabilized, the "warm water" circulation pump would be 
Started and the rate of flow to the overflow reservoir adjusted, 

4, The elevation of the overflow reservoir would be adjusted by 
using the hand winch to give tne desired sressure head as 


determined by the capacitance differential transmitter. If 
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necessary the flow rate from the "warm water" circulation pump 
would be again adjusted, 

The turbine flowmeter pulse counter and chart recorder would 

be switched on, as would the multipoint temnerature recorder, 
The circulation pump would be switched off while the ultrasonic 
unit was calibrated by means of the gauge blocks, 

The X-¥ recorder would be switched on. Tne traversing unit's 
"cycle time adjusted and noted, 

The ultrasonic transducer would be installed in the travers ing 
unit, 

All pressure taps to the manometer would de purged of air 
bubbles, The camera unit which was to record changes in the 
pressure distribution would be checked for focus and duration 
between pictures set, 

Tne coolant temperature would be set, 

The "warm water" circulation pump would be restarted, 

After the flow and temperature had again stabilized, the 
coolant pump would be turned on, 

The surface of the copper plate would be watched carefully, At 
the first appearance of an ice layer, t = t = 0, event markers 
would be manually tripped, and simultaneously the traversing 
unit and camera started, 

A test would be terminated when either the ice layer approached 
complete flow blockage, the ice layer reached a steady-state 
condition, tne motorized three way valve could no longer 
control the coolant inlet temperature, or the "warm water" 


reservoir was empted, 
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CHAPTER V 
DISCUSSION 

The study presented in this thesis is part of a long ranged investi- 
gation concerned with learning the effects caused by pnase change in 
various flow configurations and flow conditions. Until very recently 
this field has received but scant attention, Eventually it is hoped 
that enough information will be obtained to enable the development of 
a flow control device to be designed based on the freezing process, 

With this in mind the apparatus was designed to test more than the 
solution of the simplified mathematical model presented in Chapter II, 
The apparatus has been built to tilt so that the effects of natural 
convection could be looked into, The positive displacement “warm water" 
circulation pump could very simply be connected directly to the test 
section so that constant mass flow conditions could be studied. If the 
test section were to be inclined at or past the horizontal and the “warm 
water" flow rate controlled a three phase system; ice, water and air, 
could be studied-quantitatively at least. And of course the test section 
could be readily substituted for sections of various geometric configu- 
rations, 

A few worde Should be mentioned about the preliminary testing of the 
System. Tests were run at various conditions to determine the type and 
quality of the ice. Since the accuracy of the ultrasonic measurement 
technique assumes that the "test" ice be similar to the ice grown for 
Calibration purposes. WNon-homogeneities and marked variations in the 
temperature gradient across the ice modify the velocity of sound creating 


an inherent error. In general the difference between the quality of the 
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"test" and "calibration" ice varied little. Additional checks were made 
by subtracting the measured thickness of the fluid stream (time lag 
between the plexiglass wall and the ice-water interface, (see figure 
3.9) from the gap Spacing, d. Generally an accuracy of 0.005 inches 
could be expected, If, however, the ice layer grew close to the plexi- 
glass wall multiple echos would be picked up, giving senseless results, 
The cause for this seemed to be either too great a power input to the 
transducer for the close proxmity of the ice layer, or else excessive 
curvature/inclination was present in the ice-water interface, 

The ice itself grew very uniform and smooth as evidenced by the 
time lapse photographs shown in figure 5.1, Since the inlet water 
temperature possessed little superheat (to = 38°F), with a relatively 
low Reynold's number (Re =250 ), and a fairly cold coolant temperature 
(tiecl0 8 FE), thet cesinithadnitial seconde grev quite rapidly, 
eventually causing complete channel blockage. Tie ice possessing rapid 
growth is seen in the photographs as a white layer while the slower 
growing clear ice may be seen near the leading edge cf the formation, 
In general, most of the ice grown during various trial runs tended to 
be of the latter type. 

It has been intended that the design, construction, and eventual 
testing of the apparatus described in this thesis represent the first 
step in a program to research the effects caused by solidification 
occuring when the flow is through an annular passage. Fluid flow 


between parallel plates being the limiting case of such a flow, 
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APPENDIX A 
INTEGRAL CONTINUITY EQUATIGN 

Although not strictly necessary for the matheaatical formulation 
of tne problem it is worthwhile to introduce an intecral form of the 
continuity equation. Due to the change in pressure and the rate of 
solidification in the freezing zone, the mass flow is an unknown 
whicn must be determined, Thus to ensure that the mass flow is con- 
served the derived equation must be satisfied, 

Consider a control volume of unit width placed around the freezing 


zone as indicated, 


Control volume 


The mass entering the control volume between times ty and t is: 


C d 
Mi = { { PrU(Cyay a vjadyedt 
to 0 


The mass leaving the control volume in the time interval is: 


t y,(%t) 
Me = | { p U( Xe V, tidy dt 
ic 0 


6) 


> bine , enezing nt aga ans ws a 
ie i» wor? ztem sft? 13005 on re =a 


a 


“a> &b 26m auld iadt St Nae ea aunt - rin sauna rio te Ww 


- cbottetyes wd t2uft sotseapa bavtaab » ait bave a2 
fae 


wtrsogyT ames 
+s ) 


Ae Lonstg aabte stew to waite toamoo. 5 sah anoo 
-ooesto 20.9 


* 
a ¢ . A 


7 hades Syttey t 


>) 20d ooaesse Om ov foros sit puiessae ezam » a 


The fluid mass initially in the control volume is: 


Tne fluid mass in the control volume at time t is: 


X 


e 
Meg = i py, (xs t)dx 
0 


The accumulation of solid mass in the time interval is: 


Xo 
Ae 2 | R, (yO t) - y; (x; tg) )dx 
0 
Since the total mass is conserved then 


Me a Me - MM. cr Meg - Mir 


Substituting and rearranging gives: 


t d ty (x t) 
| | u(O,Y,t)dydt - | | GUX ec) Gyaat 
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Differentiating w.r.t time and formulating in terms of non dimensional 


quantities gives equation 2.3 - 9 


x n xX 
| va - [vy = 1-°? an dx 
6a ex oT 
0 0 0 


The equation introduces the additional non-dimensional parameter, p - 


tne ratio of the solid to liquid density. 


Ted fane tenentowaon fonol pints 9 pequberil aBtdsups 
ua? tot tibsen!'T o piles sii? ” ehter > 


1. ie 
ti 


APPENDIX B 


56 


57. 


- MATERIAL List 


Poor vem] nee] Sear 


a - - a _ 
= or 7 . Nt ee ee mer ee mean 


ST Te ee a te ee em meme ee ee Se ee eee ~| 
, a: ; 


‘ 


| 
; 
| 


de 


_ 


Bites! 


ss + 


ve 


= 


58 


ne ee 
BEOX — 2C/-A-O8 
a eee 
WASHES = Yyi'x 3 «Ji 
Aree i= FIEAN ©, D. * 


TT Ee a A rR 
i ‘ : 


MATERIAL LisT 


MATERIAL LIST 


TON E-= 


) ae: Pbk 5 = P 7 ¢ ¥” 5 
Se ee a ey ee " . tee. ae 
r ce ee jenn Li ae 


oa 


% sy : ms t 


"s 


MATERIAL List 


O-RING 


Ha | 
Hla) [= 
Lf Pp 


- 
g 
g 
; 


srt ig ao fe = — Sa 
ee Se See OS a 
— et 
Se ee 4 
» = how 
Cy FB 


63 


MATERIAL 


_—< 


uy \- ‘ 


See et eee 


aie 


ie 


i 6-6 


SECTIO 


ry" oS 


APPENDIX C 


64 


65 


APPENDIX C 


EQUIPMENT 


The following represents a list of the main conponents used in 


the preliminary testing of 


‘Warm Water' Storage 


Test Section 


Cooling Systen 
Coolant 

Coolant reservoir 
Refrigeration unit 


Control valve 


Coolant Pump 


Instrumentation 


Ultrasonic Test Equipment 


the system, 


150 Imptrial gallon insulated storage tank 

2 - 1500 watt Tempro Automatic water heaters 
1, ton - Copland refrigeration unit 

15 IGPH HMoyno F S Utility pump 

discharge line and return line - 1" 
diameter flexible tubing 


coolant inlet and discaarge piping - 1" 
diameter flexible tubing 

‘Warm Water’ inlet and discharge lines - 
14s" diameter flexible tubing 


Ethylene glycol-water iixture 50/50 

60 Imperial gallon capacity 

13s ton capacity - Copland 

Honeywell 1616 Three-Way mixing valve 
fiodutrol electric motor (''944) actuated 


Honeywell proportioning indicating 
temperature controller T9S4A 


Range - 50 to 100 °F 


25 IGPii Jabsco 


Branson Flow and Thickness Tester, 
Sonoray 600 
Branson Time Base Module, iiodel 610 
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Potentiometers 


Turbine Flowmeter 


Pressure Sensor 
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Honeywell Electronik 15 Strip Chart 
Multipoint Recorder 

Hewlitt Packard X-Y recorder 

Model 2 D = 2A 

Hewlitt Packard Strip Chart recorder, 
Model 7101B 


Brooks Model 4901 - 10 Journal Bearing 
design 

Brook 4300T - Frequancy Converter 
General Radio = 1192 Counter 


T, E, HM, Instrument 36 Channel 

manometer mercury filled 

Canon F-1 35 millimeter single reflex 
camera, electric motor drive and interval 
timer 

Rosemount Engineering Differential 
Pressure Transmitter Model 1151DP 
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